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Abstract — This paper considers the recovery of con- 
tinuous time signals from the magnitude of its sam- 
ples. It uses a combination of structured modulation 
and oversampling and provides sufficient conditions on 
the signal and the sampling system such that signal 
recovery is possible. In particular, it is shown that an 
average sampling rate of four times the Nyquist rate 
is sufficient to reconstruct a signal from its magnitude 
measurements. 

Index Terms — Bernstein spaces, Paley-Wiener 
spaces, phase retrieval, sampling 

I. Introduction 

In many applications, only intensity measurements are 
available to reconstruct a desired signal x. This is widely 
known as the phase retrieval problem which e.g. occurs in 
diffraction imaging applications such as X-ray crystallog- 
raphy, astronomical imaging, or in speech processing. 

In the past, concrete efforts have been made on the 
recovery of finite dimensional signals from the modulus 
of their Fourier transform. In general however, they re- 
quire strong limitations on the signal such as constraints 
■ on its z-transforms [T] or knowledge of its support [5] ■ 
Furthermore, analytic frame-theoretic approaches were 
considered in [3J, [3] which require that the number of 
measurements grows proportionally with the square of the 
space dimension. Ideas of sparse signal representation and 
convex optimization where applied in [5], [B] to allow for 
lower computational complexity. 

Note that all the above approaches addressed finite 
dimensional signals and the question is whether similar 
results can be obtained for continuous signals in infinite 
dimensional spaces. In [7] it was shown that real val- 
ued bandlimited signals are completely determined by its 
magnitude samples taken at twice the Nyquist rate. For 
complex valued signals, the results for finite dimensional 
spaces |3J, [3], [1] indicate that oversampling alone may 
not be sufficient for signal reconstruction from magnitude 
samples, since their particular choice of measurement vec- 
tors was the key to enabling signal recovery. 

In this work we are looking at continuous signals in 
Paley-Wiener spaces. Our approach extends ideas from 
E2, 0]) E] by applying a bank of modulators before the 
intensity measurements. Then we are able to reconstruct 
the signals from samples taken at a rate of four times the 
Nyquist rate. 



Basic notations for sampling and reconstruction in 
Paley-Wiener spaces are recaptured in Sec. |ITJ Sec. IIIII 
describes our sampling setup. In Sec. IIVI we provide suf- 
ficient conditions for perfect signal reconstruction from 
magnitude measurements of the Fourier Transform. The 
paper closes with a short discussion in SecfVl 



II. Sampling in Paley-Wiener Spaces 
Let § C R be an arbitrary subset of the real axis 



For 



1 < p < oo we write £ P (S) for the usual Lebesgue space 
on §. In particular, £ 2 (§) is the Hilbert space of square 
integrable functions on S with the inner product 



where the bar denotes the complex conjugate. In finite 
dimensional spaces {x, y) — y*x where * denotes the con- 
jugate transpose. Let T > be a real number. Throughout 
this paper T = [— T/2, T/2] stands for the closed interval of 
length T, and VWt/2 denotes the Paley-Wiener space of 
entire functions of exponential type T/2 whose restriction 
to R belongs to £ 2 (M). The Paley-Wiener theorem states 
that to every x £ VWt/2 there exists an x £ £ 2 (T) such 
that 



x(z) = J T x(t) e itz dt , for all z £ C . 



(1) 



In the following we will call x the signal in the time domain 
and x the signal in the Fourier domain, since its restriction 
to the real axis is a Fourier Transform. 

If not otherwise noted, our signal space will be £ 2 (T), 
i.e. we consider signals of finite energy which are supported 
on the finite interval T. This is a natural assumption since 
signals in reality are usually finite in time. It follows from 
the Paley-Wiener theorem that to every signal x £ £ 2 (T) 
there corresponds an x £ VWt/2 given by ([T]). 

A sequence A = {X n } n ei of complex numbers is said 
to be complete interpolating for VWt/2 if an d only if 
the functions {<fi n {t) ■— e _lA " t }„gz form a Riesz basis for 
£ 2 (T) 0. Let x £ £ 2 (T) be arbitrary. Then it follows from 
(□} that 



x(X n ) = (x,<p n ) 



£ 2 (T) 



for all n £ Z 



Since {<fi n }nez is a Riesz basis for £ 2 (T) the signal x can 
be reconstructed from the samples x(A) = {x(\ n )} ne z by 

x { f ) = E„ e z ( x > <M *Pn(t) = E„ e z x ( X n) i>n(t) , (2) 
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Fig. 1. Measurement setup: In each branch, the unknown signal x 
is modulated with a different sequence p( m \ m = 1, 2, . . . , M. Sub- 
sequently, the intensities of the resulting signals y( m ) are measured 
(IM) and uniformly sampled in the frequency domain. 



where {4>n}nei, is the unique dual Riesz basis of {4> n }nez 
[9]. It is well-known that in the Fourier domain 



S(z) 



S'(X n )(z-X n ) 



with S(z) =/ A lim (l- 



j1 n 

|A„|<fl " 
A„/0 



with 5a = 1 if £ A and 5j\ = otherwise, where the 
generating function S is an entire function of exponential 
type T/2. The infinite product converges uniformly on 
compact subsets of C if A is a complete interpolating 
sequence (see [IU]). Note that in order to obtain the signal 
in the time domain by (0) , it is sufficient to find the inverse 
Fourier transform of ip n (uj),uj G R. 

Example 1: The well known Shannon sampling series is 
obtained for regular sampling with A„ = n 2 ^, n G Z. 
Then S(z) = sin(-jz) and ip n ( z ) = sinc(-j[z — n^f-]) 
where sinc(x) := sin(x)/x. This corresponds to modulated 
rectangular functions in the time domain such that @ 
becomes x(t) = X^nez x(X n ) e -1 ™^"* for all tel. 

III. Measurement Methodology 

We apply a measurement methodology which uses over- 
sampling in connection with structured modulations of the 
desired signal, inspired by the approach in ;(]| . Suppose x G 
£ 2 (T) is the signal of interest. Although the loss of phase 
information is intrinsic to the measurement procedure, it 
is possible to influence the desired signal before the actual 
measurement. More precisely in our sampling scheme in 
Fig. [TJ we assume that the signal of interest x is multiplied 
with M known modulating functions p( m \ In optics, these 
modulations may be different diffraction gratings between 
the object (the desired signal) and the measurement device 
[B]. This way we obtain a collection of M representations 
(or illuminations) z/ m ) of x. Afterwards, the modulus 
of the Fourier spectra y( m ) are measured and uniformly 
sampled with frequency spacing fj. 

Let p( m ) have the following general form 



P 



„( m ) 



(3) 



where A& and a k are complex coefficients which are 
determined subsequently. Consequently the samples in the 



771th branch are given by 

= l<X„,a( m >)| 2 
with the length K vectors 



A" 



5>i m) x(n[3 + X k ) 



k=l 



(4) 
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and x„ 



/ x{n0 + Ai) 



\ x{n[3 + X K ) 



We will show that if a' 7 ™-* and the interpolation points 
A n ,fc := n/3 + Xk : n G Z ,k = 1,...,K are properly 
chosen ; it is possible to reconstruct x from all samples 
c = {cl m) : m = 1,. ..,M; n G Z}. 

A . Choice of the coefficients ct k m ^ 

In order to determine the vector x„ G from 
the M intensity measurements c( m ), we apply a re- 
sult from [Ij. It states that if the family of C^-vectors 
A = {a^, . . . , a( M '} constitutes a 2-uniform Af/AT-tight 
frame which contains M = K 2 vectors or A is a union of 
K+ 1 mutually unbiased bases in C , then every x n G C 
can be reconstructed up to a constant phase from the 
magnitude of the inner products ([?]). For simplicity, we will 
only discuss the first case here and therefore fix M = K 2 . 

Condition A: A sampling system as in Fig. [7] is said to 
satisfy Condition A if the coefficients ce k " 1 ^ in Q are such 
that A constitutes a 2-uniform M / K -tight frame. 

Then reconstruction will be based on the following formula 
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with rank-1 matrices Q x = xx*. For K = 2 a valid choice 
for A reads [J 



,(2) 



,(3) 
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with a 



1(1 - -1=) and (3 = e i5w / 4 



B. Choice of the interpolation points 

Let {Afc}^L 1 be ordered increasingly by their real parts. 
For each n G Z, the vector x„ contains the values of x 
at K distinct interpolation points in the complex plane 
collected in the sequences 

K--={KJti with X a n<k =np + \ k , neZ. (6) 

Therein, the parameter a G N denotes the number of 
overlapping points of consecutive sets (cf. also Fig|2). 
More precisely, we require for every n G Z that 



\ a \a 

A n,i ~~ A n-1, 



K-i+l 



for alii = 1, 



,a . 



(7) 
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Fig. 2. Illustration for the choice of interpolation points in the 
complex plane for K = 6 in (O and an overlap a = 2. 

In the following we denote by A n = A" n A^ +1 the set 
of overlapping interpolation points between A° and A^ +1 , 
and we define the overall interpolation sequences 

A a := U„ e zA° . 

In general we allow for a > 1, but we will see that a = 1 
is generally sufficient for reconstruction. 

As explained in Sec. [TTJ x £ £ 2 (T) can be perfectly 
reconstructed by @ if A a is complete interpolating for 
VW T/2 . 

Condition B: A sampling system as in Fig. Q] is said 
to satisfy Condition B if the coefficients {AfcjjfLj in ([3]) 
are such that A a is complete interpolating for VWt /2 an d 
satisfies ([7]) for a certain 1 < a < K . 

Note that interpolating sequences fulfilling Condition B 
are /3-periodic. In general it is hard to characterize sets 
which fulfill this condition. However one specific example 
is when A a equals to the set of zeros of a /3-periodic sine- 
type function of type T'/2 < T/2 (see, e.g., [§], [TT]). Sine- 
type functions are entire functions / of exponential type 
T'/2 with simple and isolated zeros and for which there 
exist positive constants A, B, H such that 

Ae% M < + i v )\ < Be^ 1 , for |n| > H . 

Moreover, sequences which differ from interpolation pat- 
terns derived from sine-type functions only by their imag- 
inary parts, are still zeros of sine-type functions [12] . The 
complete interpolating property is also preserved under 
little shifts in the real part using Katsnelson's theorem 
[llj . Apparently it is possible to construct many non- 
uniform complex interpolation sequences A a which satisfy 
Condition B. One particular simple construction is ob- 
tained by starting with the zeros of the sine-type function 
sin(^-z), which has equally spaced zeros on the real axis 
(cf. Example [JJ. 

IV. Signal Reconstruction - Main Results 

We assume a sampling scheme as described in Section HIIl 
(cf. Fig. [T]) which satisfies Condition A and B. For this 
setup, we show that basically every x £ £ 2 (T) can be 
reconstructed from the samples in ^ . The proof provides 
an explicit algorithm for perfect signal recovery. 



Theorem 1: Let x £ C 2 (T) be sampled according to the 
scheme in Section \IIH which satisfies Condition A and B, 
and let c = {cl m ''}™^ 1 '"'' M be the sampling sequence in 
If the set x{A n ) contains at least one non-zero element 
for each n £ Z, then x can be perfectly reconstructed from 
c up to a constant phase. 

Proof: According to Condition B of the sampling 
system, A a is complete interpolating for VWt/2- There- 
fore the signal x can be reconstructed from the vectors 
{x n }„gz using It remains to show that {x„} n gz can 
be determined from c. 

Let n £ Z be arbitrary. Since the sampling system 
satisfies Condition A, we can use (0) to obtain the rank-1 
matrix Q n := x„x* from the measurements {c£ }^f =1 . 
Then x„ £ C K is obtained by factorizing Q n . However, 
such a factorization is only unique up to a constant phase 
factor. If the phase <f> n ^ of one element [x„]i is known, the 
vector x„ can be completely determined from Q n by 

x(n[3 + Afe) = > /[Qn]Me I( *"' i " Wg(l ° B]l '* )) . Vfc ± 1 ■ ( 8 ) 

Assume that we start the determination of the sequence 
{x n } Iig z at a certain hq £ Z. In this initial step, we set 
the constant phase of x„ arbitrarily. In the next step, 
we determine x„ 0+ i. After the factorization of <3n +i, the 
vector x, l0 +i is only determined up to a constant phase. 
However, since Aq ng is non-empty, and because x(Aq ) 
contains at least one non-zero element, we have phase 
knowledge of at least one entry of x„ 0+ i, say a;(A° +1 A ), 
where A^ 0+1 i is an overlapping interpolation point of A^ o 
and A° o+1 . Thus, we can completely determine x„ 0+ i and 
successively all n = no ± 1, no ± 2, . . . using (JSJ to obtain 
x(A a )e i9 °. 

Note that the arbitrary setting of the phase of the initial 
vector x„ yields a constant phase shift 8q £ [— ir, it) for 
all x„ which persists also after the reconstruction of the 
time signal by |2}. ■ 

Theorem Q] states that x £ C 2 (T) can only be re- 
constructed if x £ VWt/2 has at most a — 1 zeros on 
the overlapping interpolation sets Aq n . However, this 
restriction is not too limiting, since it is known that the 
zeros of an entire function of exponential type can not be 
arbitrarily dense in C. For example, defining Z n :— {z £ 
C : rnr/T < \z\ < (n+l)-ir/T}, the result in [13] states that 
for every x £ VWt/2 there exist only finitely many sets 
Z n which contain more than one zero of x. Consequently, 
by choosing the spacing of the interpolation points in the 
overlapping sets Aq less than tt/T, it is very unlikely that 
a randomly chosen function from 'PVVt/2 fails to satisfy 
the condition of Theorem [TJ especially for a > 1. 

To avoid such pathological cases, we may a priori restrict 
the function space allowed in Theorem [T] to prevent x 
from having zeros in A a . To this end, we first state an 
extension of a well-known lemma by Duffin, Schaeffer }14j . 
The straightforward proof is omitted here. 



Lemma 2: Let x(z) G VWt/2 be an entire function of 
z = £ + in satisfying \x(^)\ < M on the real axis. Then for 
every T' > T the function 

y{z) = Mcos(yz) - x{z) (9) 

belongs to the Bernstein space B^, ,„ and there exists a 
constant H = H{T,T') such that \y(z)\ > Vz : \rj\ > H. 

The Bernstein space B%?,/ 2 * s the set of all entire functions 
of exponential type T"/2 whose restriction to R is in 
£°° (R) . Upon this we can establish the following corollary 
for functions in K M (T) = {x G C 2 (T) : ||ac|Ui( T ) < M}. 
Corollary 3: Let x G TZm{^) be sampled according to the 
scheme in Sec. UJil Then there exist interpolation sequences 
A a with overlap a > 1 such that every x G 7£m(T) can be 
perfectly reconstructed (up to a constant phase) from the 
measurements (HJ). 

Sketch of proof : The bounded £ 1 -norm of x implies that 
|Zc(f)| < M for all £ G K. Fix T' > T, then y, de- 
fined by (O, has zeros only for 1 77 < H by Lemma [2] 
Choose A a as the zero set of a sine-type function of 
type T'/2. Then by [T2] we can shift the imaginary parts 
of the interpolation points such that > H for all 
k while A a remains to be the zero set of a sine-type 
function. Now assume that the measurement is preceded 
by an addition of a cosine as in © such that as in 
the proof of Theorem [TJ the factorization step yields 
y„ := [y(n[3 + Ai), . . . ,y(nf3 + A K )] T = Mu n - x„ with 
u„ := [cos(f [n/3 + Ai]), . . . ,cos(~ [nfi + Xk])] t ■ We then 
obtain y(A°)e ie ° since the overlapping interpolation points 
do not coincide with zeros of x and the phase information 
can be propagated. Since y G &t'/2 an ^ because A a is 
the zero set of a sine-type function, we know by [15, The- 
orem 2] that it is possible to stably reconstruct y(z) e 
from y(A a ) e w ° using y(z) = J2 n ezV( X n)^n(z). Note that 
it is not possible to simply extract x using ^ since #0 is 
unknown. However we can determine 

x'{z) = M cos(^z) - y(z) e i9 ° 

= x(z)e i9 ° +M(1- e i9o )cos(^z). 

Applying the inverse Fourier transform to x! yields x' (t) = 
x(t)e w ° yt G T, using that T > T. It is therefore possible 
to reconstruct the time signal x up to a constant phase. ■ 

V. Discussion and Outlook 

To determine the sampling system in Fig[TJ one has to 
fix K, M, a and j3. The number K > 2 can be chosen 
arbitrarily. Then M — K 2 is fixed, and 1 < a < K — 1. The 
sampling period f3 has to be chosen such that the sampling 
system satisfies Condition B and in particular that A a is 
complete interpolating for VWt/2- As discussed before, 
one possible choice may start with the zeros of the function 
sin(^z) with T > T. Then 8 := \ k - A fe _i = 2tt/T' such 
that (3 = (K — a)S. Therewith, the total sampling rate 



becomes 

, M K 2 K 2 T K 2 V 
R{a,K) = — = — — = — — = — — R Ny 

p (K — a)o K — ciztt K — a 1 

where R^ y := T/(2tt) is the Nyquist rate. It is apparent 
that R(a, K) grows asymptotically proportional with K 
and the oversampling factor T'/T and increases with the 
overlap a. The minimal sampling rate R(a, K) = 4i?N y 
is obtained for K = 2, T = T, and a = 1. Since T'/T 
can be made arbitrarily close to 1, and in connection with 
Theorem [T] and Cor. [3] this shows that a sampling rate 
greater 4i?Ny is sufficient for signal recovery based on the 
magnitude of the signal samples. This corresponds to the 
findings in [3] for finite dimensional spaces, where it was 
shown that basically any x G C N can be reconstructed 
from M > AN — 2 magnitude samples. 

Finally, we note that the above framework can be 
applied exactly the same way for bandlimited signals. To 
this end, one only has to exchange the time and frequency 
domain. Then the modulators in Fig.[T]have to be replaced 
by linear filters and the sampling of the magnitudes has 
to be done in the time domain. 
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